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We investigate the charge- and spin dynamical structure factors for the 2D one-band Hubbard 
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doping, our results provide a consistent overall picture of the interplay between charge, spin and 
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I. INTRODUCTION 

Two-particle correlation functions, such as the dynam- 
ical spin- and charge correlation functions, determine a 
variety of crucial properties of many-body systems. Their 
poles as a function of frequency and momentum describe 
the elementary excitations, i.e. electron-hole excitations 
and collective modes, such as spin- and charge-density 
waves. Furthermore, an effective way to identify contin- 
uous phase transitions is to search for divergencies of sus- 
ceptibilities, i.e. two-particle correlation functions. Yet, 
compared to studies of single-particle Green's functions 
and their spectral properties, where a good overall ac- 
cord between theoretical models (Hubbard type-models) 
and experiment (ARPES) has been established (see Refs. 

the situation is usually not so satisfying for two- 
particle Green's functions. This is especially so for the 
case of correlated electron systems such as high- Tc su- 
perconductors (HTSC). The primary reason for this is 
that calculations of these Green's functions are, from a 
numerical point of view, much more involved. 

To expose the problem let us consider the spin- 
susceptibility which is given by: 



^(2) = Jl^^hA^) .with 
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Here, L corresponds to the lattice size, /3 is the inverse 
temperature and q = (gr,l]^), q being the momentum 
and ftrn a (bosonic) Matsubara frequency. To simplify 
the notation, we have adopted a path integral coherent 
state notation with Grassman variables: 
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Bethe-Salpeter equation for the two-particle propa- 



The two-particle irreducible vertex, F^/ /^//(g), is defined 



through the Bethe-Salpeter equation, 

k\k" 

which is diagrammatically depicted in Fig. [H 

Within the Dynamical Cluster Approximation 
(DCA)^i^, which we consider in this paper, one can 
consistently define the two-particle Green's functions, 
by extracting the irreducible vertex function from the 
cluster. This approximation maps the original lattice 
problem to a cluster of size Lc = Ic x Ic embedded 
in a self-consistent host. Thus, correlations up to a 
range < Ic are treated accurately, while the physics at 
longer length-scales is described at the mean-field level. 
The DCA is conveniently formulated in momentum 
space, i.e. via a patching of the BZ. Let K denote 
such a patch, and k the original lattice momentum. 
The approximation boils down to restricting momentum 
conservation only to the patch wave vectors K. This 
approximation is justified if /c-space correlation func- 
tions are rather structureless and, thus, in real space 
short-ranged. From the technical point of view, the 
approximations are implemented via the Laue function 
A(/ci, /c2, ^3, ^4), which guarantees momentum conser- 
vation up to a reciprocal lattice vector. In the DCA, the 
Laue function is replaced by A]jca{Ki^ Ks^ 
thereby insuring momentum conservation only between 
the cluster momenta. It is understood that ki belongs 
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to the patch Ki. 

To define uniquely the DCA approximation, in par- 
ticular in view of two-particle quantities, it is useful to 
start with the Luttinger-Ward functional <l>, which is 
computed using the DCA Laue function. Hence, ^dca 
is a functional of a coarse-grained Green's function, 
G{K ^iujrn) = G(K). Irreducible quantities such as the 
self-energy, and the two-particle irreducible vertex are 
calculated on the cluster and correspond, respectively, 
to the first- and second-order functional derivatives of 
^DCA with respect to G. Using the cluster irreducible 
self-energy, S(i^), and two-particle vertex, Tx' ,k"{Q)^ 
one can then compute the lattice single-particle and lat- 
tice two-particle correlation functions using the Dyson 
and Bethe-Salpeter equations. This construction of two- 
particle quantities has the appealing property that they 
are thermodynamically consistent^'^. Hence, the spin 
susceptibility as calculated using the particle-hole cor- 
relation functions corresponds precisely to the derivative 
of the magnetization with respect to an applied uniform 
static magnetic field. The technical aspects of the above 
program are readily carried out for single-particle prop- 
erties. However a full calculation of the irreducible two- 
particle vertex — even within the DCA — is prohibitively 
expensive^ and, thus, has never been carried out. 

In the present work, we would like to overcome this sit- 
uation by suggesting a scheme where the K^— and K" — 
dependencies of the irreducible vertex are neglected. At 
low temperatures, this amounts to the assumption that 
in a energy and momentum window around the Fermi 
surface, the irreducible vertex depends weakly on K^ — 
and Following this assumption, an effective two- 

particle vertex in terms of an average over the K^— and 
K^^ dependencies of rx',K"(Q) is introduced: 



1 

PL 



Ueff{Q) = {TK',K<Q))- 



(4) 



As shown in an earlier Quantum Monte Carlo (QMC) 
study by Bulut et al. (Ref. ITol ) for a single QMC cluster, 
this is reasonable for the 2D Hubbard model (on this 
QMC cluster of size 8x8 with U = 8t). The authors of 
Ref. dO have also calculated with Ueff{Q) the effective 
electron-electron interaction for the 8x8 single QMC 
cluster. Both the momentum and frequency dependence 
were in rather good agreement with the QMC results for 
the effective electron-electron interaction. 

Replacing the irreducible vertex by -^Ueff{Q) in the 
cluster version of the Bethe-Salpeter Eq. ([3|) and carrying 
out the summations to obtain the cluster susceptibility 
gives: 



Ueff iO) = 
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1 



xo(Q) x(Q)' 



(5) 



where x corresponds to the fully interacting susceptibility 
on the DCA cluster in the particle-hole channel and xo 
is the corresponding bubble as obtained from the coarse- 
grained Green's functions. Finally, our estimate of the 



lattice susceptibility reads: 
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where Xo{q) corresponds to the bubble of the dressed 
lattice Green's functions: G{k) = i s-ri^^ - 

In the following section, we describe some aspects of 
the explicit implementation of the DCA which is based 
on the Hirsch-Fye QMC algorithm^^. Our new approach 
requires substantial testing. In Sec. IIII Al we compare the 
Neel temperature as obtained within the DCA without 
any further approximations to the result obtained with 
our new approach. In Sec. IIIIB| we present results for 
the temperature and doping dependence of the spin and 
charge dynamical structure factors and compare the high 
temperature data with auxiliary field QMC simulations. 
Finally, Sec. IIVI draws conclusions. 



II. IMPLEMENTATION OF THE DCA 

We consider the standard model of strongly correlated 
electron systems, the single-band Hubbard mode\^^^ 

H = -^Ujclcj^ + U ^m^nn, (7) 



with hopping between nearest neighbors t and Hubbard 
interaction U. The energy scale is set by t and through- 
out the paper we consider U = St. 

Our goal is to compute the spin, S{q^uj) and charge 
C{q^uo) dynamical structure factors. They are given, re- 
spectively, by: 

{S\q, T)S'i-q, 0)) = ll dw e"-'^ S(q, u) (8a) 

{N{q, T)N{-q, Q)) = \j dw e"™ C{q, to). (8b) 

Here, S^{q) = (n,-,T - n,a) and N{q) = 

Xlj 6*^^ (^j,T + ' "^^^ ^^^^ hand side of the above 
equations are obtained from the corresponding suscepti- 
bility as calculated from Eq. (|6]). Finally, a stochastic 
version of the Maximum Entropy method^^i^ is used to 
extract the dynamical quantities. 

In order to cross check our results, we slightly modify 
Eq. 16]) to: 
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Here, we have introduced an additional "controlling" pa- 
rameter a in the susceptibility denominator, which is cal- 
culated in a self-consistent manner. It assures, for exam- 
ple in the case of the longitudinal spin response, that x(g) 
obeys the following sum rule (a similar idea, to use sum- 
rules for constructing a controlled local approximation for 
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FIG. 2: Static {i^m = 0) irreducible particle- hole interaction 
Ueff for different cluster momentum vectors and dopings at 
inverse temperature pt = 6 and U = St. 



the irreducible two-particle vertex has been implemented 
by Vilk and Tremblay (Ref. [13)): 
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Of course, a should be as close as possible equal to a = 1, 
which is indeed what we will find after implementing the 
sum rule (see below). 

Our implementation of the DCA for the Hubbard 
model is standard^. Here, we will only discuss our inter- 
polation scheme as well as the implementation of a SU(2)- 
spin symmetry broken algorithm. Since the DCA evalu- 
ates the irreducible quantities, ^{K_) as well as Ueff{Q) 
for the cluster wave vectors, an interpolation scheme has 
to be used. To this, we adopt the following strategy: for 
a fixed Matsubara frequency iQrn and for each cluster 
vector Q, the effective interaction Ueff is rewritten as a 
series expansion: 

UeffiQ^inm) = ^^e^^^^A,(il]^), (11) 

i Ai 

with i = 0,...,n — 1, where n is the number of the clus- 
ter momentum vectors Q. The quantity represents 
vectors, where each vector from the corresponding 
belongs to the same "shell" around the origin (0,0) in 
real space, i.e. 
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With a given Uef / , Eq. (pT]) can be inverted to uniquely 
determine Ai. With these coefficients, one can compute 
the effective particle-hole interaction for every lattice mo- 
mentum vector q. This interpolation method works well 
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FIG. 3: SU(2) symmetry broken DCA calculation. Left: A.F. 
unit cell with new basis vectors in realspace. Right: Reduced 
A.F. unit cell in momentum space. 



when Ueff is localized in real space and the sum in Eq. 
(pTj) can be cut-off at a given shell. 

The effective particle- hole interaction Uef / is shown in 
Fig. [2] for a variety of dopings at inverse temperature 
Pt = 6^ U/t = 8 and on an Lc = 8 cluster, which cor- 
responds to the so-called "8A" Betts cluster (see Refs. 
[isIIiqI ). The /7e// -function displays a smooth momen- 
tum dependence. These observations further support the 
interpolation scheme (Eq. [11]). Thus, indeed, Ueff is 
rather localized in real space with sizable reduction from 
its bare U = St value for larger doping and a further 
slight reduction at gr = (tt, tt). The reduction is partly 
due to the self-energy effects in the single-particle propa- 
gator, which reduce xo from its non-interacting {U = 0) 
value y^^^ . Partly, it also reflects both the Kanamori (see 
Ref. UM repeated particle-particle scattering and vertex 
corrections. 

Summarizing, the new approach to two-particle prop- 
erties relies on two approximations which render the cal- 
culation of the corresponding Green's function possible. 
Firstly, the effective particle-hole interaction Ueff{Q) de- 
pends only on the center-of-mass momentum and fre- 
quency, i.e. Q and iVtrn- Secondly, x{Q)i is extracted 
directly from the cluster and Xo(Q) is obtained from the 
bubble of the coarse-grained Green's functions. 

To generate DCA results for the Neel temperature, we 
have used an SU(2) symmetry broken code. The setup is 
illustrated in Fig. [3l We introduce a doubling of the unit 
cell — to accommodate AF ordering — which in turn 
defines the magnetic Brillouin zone. The DCA k-space 
patching is carried out in the magnetic Brillouin zone and 
the Dyson equation for the single-particle propagator is 
given as a matrix equation: 



with 



G%k) = 



(12) 



(13) 



With the SU(2) symmetry broken algorithm, one can 
compute directly the staggered magnetization, i.e. m = 
e^^^ {rij^^ — Tij^l)^ and thereby determine the tran- 
sition temperature. Since the DCA is a conserving ap- 
proximation, the so determined transition temperature 



4 



0.6 
0.5 



0.4 - 
0.3 
0.2 
0.1 




vv 

V V 



O o 
vv V 



A 
A 



A 
A 
A 



n \ 

PM-phase (symm. breaking) Q 

AF-phase (symm. breaking) 

PM-phase (2-particle) ^ 

AF-phase (2-particle) A 
indeterminable 



O V VV V V 

O lDD V V V T- 

^ VVV V 
A A O O A V ▼ 
A 0. A 

_J \ 



0.05 



0.1 

5 



0.15 



0.2 



FIG. 4: Phase diagram for the one-band Hubbard model with 
U = 8t with respect to different temperatures and fiUings. 
The calculations are carried out on an Lc = 8 cluster. The 
red and gray (blue and blank) objects indicate the antifer- 
romagnetic (paramagnetic) phase. Shading of AF and PM 
regions is a guide to the eye. Details are in the text. 
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FIG. 5: Static irreducible particle-hole interaction Ueff for 
the cluster momentum vector Q = (tt, tt). The inset shows the 
static free lattice susceptibility xo for the momentum vector 
q — (tTjTt). The bare Hubbard interaction strength \sU — 8t. 



corresponds precisely to the temperature scale at which 
the corresponding susceptibility, calculated without any 
approximations on the irreducible vertex Tx'^x"{Q)^ di- 
verges. 



III. RESULTS 

A. Comparison of the AF transition temperature 
with a symmetry broken DCA calculation 

A first test of the validity of our new approach is a 
comparison with the SU(2) symmetry broken DCA cal- 
culation on an Lc = 8 cluster at t/ = St. The idea is 
to extract the Neel temperature Tn from a divergence 
in the spin susceptibility as calculated in the above de- 
scribed (paramagnetic) scheme — see Eq. (|9]) — and to 
compare it to the DCA result as obtained from the SU(2) 
symmetry broken algorithm. This comparison provides 
information on the accuracy of our approximation to the 
two-particle irreducible vertex (see Eq. [4j). 

Using the SU(2) symmetry breaking algorithm, the 
magnetic phase diagram for the one-band Hubbard model 
as a function of doping is shown in Fig. [H The para- 
(antiferro) magnetic phase transition is indicated here by 
gray (blank) circles. At half-filling T/v — OAt and mag- 
netism survives up to approximately 15 % hole doping. 
It is know that the convergence of the magnetization dur- 
ing the self-consistent steps in the DCA approach is ex- 
tremely poor near the phase transition and, therefore, 
we cannot estimate the transition temperature more pre- 
cisely than shown in Fig. [H However, our precision is 
sufficient for comparison. We again stress that the so 
determined magnetic phase diagram corresponds to the 
exact DCA result where no approximation — apart from 
coarse graining — is made on the particle- hole irreducible 



vertex. 

The blue (red) triangles indicate the transition line for 
the para- to the antiferromagnetic solutions extracted 
from the divergent spin susceptibility (Eq. [9]) within 
the paramagnetic calculation. A precise estimation of 
the Neel temperature requires very accurate results and 
boils down to finding the zeros of the denominator of 
Eq. ([9]). In Fig. O we consider the effective irreducible 
particle- hole interaction Ueff for the static case and for 
the cluster momentum Q = (tt, tt) relevant for the AF 
instability. As apparent, the irreducible particle- hole in- 
teraction becomes weaker with increasing doping. On 
the other hand, the susceptibility xo(^7^^m = 0) grows 
with increasing doping. At a first glance both quantities 
Ueff and Xo (see Fig. [5j) vary smoothly as a function of 
doping. However, in the vicinity of the phase transition, 
signalized by the vanishing of the denominator in Eq. (|9]), 
the precise interplay between Uef / and xo becomes deli- 
cately important and renders an accurate estimate of the 
Neel temperature difficult. Given the difficulty in deter- 
mining precisely the Neel temperature, we obtain good 
agreement between both methods at (5 > 10 %. Note 
that in those calculations the values of a 0.86 — 0.97 
are required to satisfy the sum rule in Eq. (p!Q|) . At 
smaller dopings, and in particular at half-band filling the 
Neel temperature, as determined by the vanishing of the 
denominator in Eq. (|9]), underestimates the DCA result. 
Hence, in this limit, the and K^' dependence of the 
irreducible vertex plays an important role in the deter- 
mination of Tn and cannot be neglected. 

Let us emphasize, that a good agreement between the 
Neel temperatures at J :^ 10 % and above is a non-trivial 
achievement lending substantial support to the above 
new scheme for extracting two-particle Green's functions. 
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FIG. 6: DCA (left) versus auxiliary field QMC (BSS) (right) 
for the dynamical spin and charge structure factors of the 
Hubbard model at U/t = 8, 6 ^ U % and pt = 3. The 
BSS data on the 8x8 lattice is essentially exact and acts as 
a benchmark for the DCA approach. The DCA calculations 
were carried out on an Lc = 8 cluster. Here we have used 
a = 0.98 and a = 1.01 to satisfy the sum rule in the spin and 
charge sectors respectively. 
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FIG. 7: DCA a) versus BSS b) static spin correlation function 
ai U / 1 = 8, 6 ^ 14 % and /3t = 3 on an Lc = 8 cluster. 



B. Dynamical Spin and Charge structure factors 

To further asses the validity of our approach, we com- 
pare it to exact auxiliary- field Blankenbecler, Scalapino, 
Sugar (BSS) QMC results (Ref. i). This method has a 
severe sign-problem especially in the vicinity of (5 :^ 10 % 
and, hence, is restricted to high temperatures. 

Such a comparison for the dynamical spin, S{q^ cj), and 
charge, C{q^uo) dynamical structure factors is shown in 
Fig. Elat = 3, (5 14 % and U/t = 8. The BSS results 
correspond to simulations on an 8 x 8 lattice. Fig. [Bb) 
depicts the BSS-QMC data in the spin sector. Due to 
short-range spin-spin correlations, remnants of the spin- 
density- wave are observable, displaying a characteristic 
energy-scale of 2 J, where J is the usual exchange cou- 
pling, i.e. J = Ajj. The two-particle DCA calculations 
show spin excitations with the dominant weight concen- 
trated, as expected and seen in the QMC data, around 
the AF wavevector (7r,7r). As apparent from the sum- 
rule, 

{S%q)S%-q)) = ^ ldu;Siq,L0) (14) 

(see Fig. [7] a)), the DCA overestimates the weight at this 
wave vector but does very well away from q = (7r,7r). 
The dispersion in the two-particle data has again a higher 
energy branch around 2 J, but it also shows features at J. 
Since the total spin is a conserved quantity, one expects 
a zero-energy excitation at gr = (0,0). This is exactly 



reproduced in the 8x8 QMC-BSS data, and qualitatively 
in the DCA results. 

As a function of decreasing temperature, the DCA dy- 
namical spin structure factor shows a more pronounced 
spin- wave spectrum. This is confirmed in Fig. [8] on the 
left hand side. Here, we fix the temperature to /3t = 6 
and keep the doping at (5 14 % but vary the clus- 
ter size. As apparent, for all considered cluster sizes 
{Lc = 4, 8, 10, 16) a spin wave feature is indeed observ- 
able: a peak maximum at g = (tt, tt) is present and the 
correct energy scale at gr = (tt, 0) with 2 J is recovered. 
Unfortunately, a direct comparison of both calculations 
at lower temperature is not possible due to the severe 
minus-sign problem in the BSS calculation. 

The investigation of the dynamical charge correlation 
function for the above parameters shows that the DCA 
calculations, which are depicted in Fig. [6] c), can also 
reproduce basic characteristics of the BSS charge excita- 
tion spectrum [6] d) . Both calculations show excitations 
at ~ which are set by the remnants of the Mott- 
Hubbard gap. Similar results are obtained at lower tem- 
peratures {l3t = 6) on the right hand side of Fig. [8] for 
different cluster sizes {Lc = 4, 8, 10, 16). The correspond- 



TABLE I: Values of a for the spectra in Fig. [H 
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FIG. 8: Dynamical spin and charge structure factors of the 
Hubbard model at (3t = 6, 6 ^ 14 % and U/t = 8. for different 
cluster sizes: (a-b): Lc = 4, (c-d): Lc — 8, (e-f): Lc = 10 and 
(g-h): Lc = 16. 



ing values of a are listed in Tab. [H These values confirm 
the overall correctness of our approach in that the corre- 
sponding sum rule for the charge response is accurately 
(exactly for a = 1) fulfilled. 

The doping dependence of the spin- and charge- 
response is examined in Fig [9l Here, we restrict our 
calculations to the Lc = S cluster at /3t = 6 and dop- 
ings between ^ = 14 % and S = 32 %. At ^ = 14 % 
(see Fig. [8]) the dynamical spin structure factor dis- 
plays a spin wave dispersion with energy scale J. That 
is E^^^{7T, 0) = 2 J with J = 4|^. As the system is fur- 
ther doped {S = 27 %) the dispersion is no longer sharply 
peaked around q = (tt, tt). The excitations broaden up 
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FIG. 9: Spin- and charge structure functions for different dop- 
ings: (a,b): 27 % and (c,d): 32 % The calculations are carried 
out on an Lc = 8 cluster at (3t = 6. 



and change their energy scale from J = 4|j to an en- 
ergy scale set by the non-interacting bandwidth. This 
effect becomes even more visible with higher dopings at 
S = 32 % (Fig. [9]c)). Furthermore, the spectrum of the 
charge response shows a reduction of the weight of states 
at high energies {uj/t ^ 8). This behavior corresponds to 
the loss of weight of the upper Hubbard band with in- 
creasing doping. The corresponding equal time spin and 
charge correlation functions of Fig. [9] c-d) are depicted 
in Fig. [TOl As in auxiliary-field QMC simulations^, the 
equal time spin correlation function shows a set of peaks 
at gr = (tt ± e, tt) and q = (tt, tt ± e). Here e is propor- 
tional to the doping. The overall trend of the doping 
dependence of the spin- and charge-responses is in good 
agreement with the previous findings of QMC simulations 
(Refs.lHjsl): there it was shown that the spin-response has 
a characteristic energy scale a; ~ 2 J and an SDW-like dis- 
persion up to about S ^ 10 — 15 % doping, despite the 
fact that at these dopings the spin-spin correlations are 
very short-ranged (of order of the lattice parameter). 

A lot of the features of the two-particle spectra have 
direct influence on the single-particle spectral function 
and vice-versa. At optimal doping, S = 14 % the spec- 
tral function A{q^ co) in Fig. [TT]a) shows three distinguish 
features. An upper Hubbard band {uo/t ^ St) and a lower 
Hubbard band which splits in an incoherent background 
and a quasiparticle band of width set by the magnetic 
scale J. In agreement with earlier QMC data (Refs.[2|0), 
we view this narrow quasiparticle band as a fingerprint of 
a spin polaron where the bare particle is dressed by spin 
fluctuations. The fact that the dynamical spin structure 
factor in Fig [8] c) shows a well defined magnon desper- 
sion at this temperature and doping, (5 = 14 %, allows 
us to interpret the features centered around q = (0,0) 
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FIG. 10: Static spin a) and charge b) correlation function at 
U/t = 8, ^ 32 % and /3t = 6 on an Lc = 8 cluster. 



and below the Fermi energy as backfolding or shadows 
of the quasiparticle band at = (tt, tt). A comparison of 
the charge response spectrum in Fig. [8] d) with the cor- 
responding single-particle spectra in Fig. [TT] a) reveals 
that the response in the particle-hole channel at almost 
zero energy is caused by particle-hole excitations around 
the quasi-particle spin-polaron band close to the Fermi 
energy. The high energy excitations, mentioned above, 
are due to transitions from the quasi-particle band to the 
upper Hubbard band. 

As a function of doping, notable changes in the spectral 
function which are reflected in the two-particle proper- 
ties are apparent. On one hand, the spectral weight in 
the upper Hubbard band is reduced. As mentioned pre- 
viously, this reduction in high energy spectral weight is 
apparent in the dynamical charge structure factor. On 
the other hand, at higher dopings the magnetic fluctu- 
ations are suppressed. Consequently, the narrow band 
changes its bandwidth from the magnetic exchange en- 
ergy J to the free bandwidth. This evolution is clearly 
apparent in Figs. [TT]b) and c) and is in good agreement 
with previous BSS-QMC results^. 



IV. CONCLUSIONS 

Two-particle spectral functions, such as spin- and 
charge- dynamical spin structure factors, are clearly 
quantities which are of crucial relevance for the under- 
standing of correlated materials. Within the DCA, those 
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FIG. 11: Angle-resolved spectral functions A{q^ijj) for various 
hole dopings: a): 14 %, b): 27 % and c): 32 %. Calculations 
are carried out on an Lc = 8 cluster at f3t — 



quantities have remained elusive due to numerical com- 
plexity. In principle, the two-particle irreducible vertex, 
^ K' ,K"{Q)i containing three momenta and frequencies 
has to be extracted from the cluster and inserted in the 
Bethe-Salpeter equation. Calculating this quantity on 
the cluster is in principle possible, but corresponds to a 
daunting task which — to the best of our knowledge — 
has never been carried out. To circumvent this problem, 
we have proposed a simplification which relies on the 
assumption that rx',K"(Q) is only weakly dependent on 
and K^' . Given the validity of this assumption, one 
can average over and K^' and retain its dependency 
on frequency and momenta, Q, of the excitation. We 
have tested this idea extensively for the two dimensional 
Hubbard model at strong couphngs, U/t = 8. At dopings 
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6 ^ 10 % and higher, we have found a good agreement 
between the Neel temperature on an Lc = 8 lattice 
as calculated within a symmetry broken DCA scheme 
with our new two-particle approach. This finding lends 
support to the validity of our scheme in this doping 
range. Furthermore, we studied the doping and tem- 
perature dependence of the spin and charge dynamical 
structure factors as well as the single-particle spectral 
function at f3t = 6. Our results provide a consistent 
picture of the physics of doped Mott insulators, very 
reminiscent of previous findings within auxiliary field 
QMC simulations. The strong point of the method, in 
contrast to auxiliary field QMC approaches, is that it 
can be pushed to lower temperatures above and below 
the superconducting transition temperature. Work in 
this direction is presently under progress. 
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